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Abstract. Dense suspensions of self-propelled rod- like particles exhibit a fascinating 
variety of non-equilibrium phenomena. By means of computer simulations of a 
minimal model for rigid self-propelled colloidal rods with variable shape we explore 
the generic diagram of emerging states over a large range of rod densities and aspect 
ratios. The dynamics is studied using a simple numerical scheme for the overdamped 
noiseless frictional dynamics of a many-body system in which steric forces are dominant 
over hydrodynamic ones. The different emergent states are identified by various 
characteristic correlation functions and suitable order parameter fields. At low density 
and aspect ratio, a disordered phase with no coherent motion precedes a highly- 
cooperative swarming state at large aspect ratio. Conversely, at high densities weakly 
anisomctric particles show a distinct jamming transition whereas slender particles form 
dynamic laning patterns. In between there is a large window corresponding to strongly 
vortical, turbulent flow. The different dynamical states should be verifiable in systems 
of swimming bacteria and artificial rod-like micro-swimmers. 
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1. Introduction 

Collections of swimming microorganisms and self-propelled particles are able to form 
remarkable macroscopic patterns pfl El [3j H] including swarms [3 [6J and complex 
vortices El El EH E] • The tendency for neighbouring particles to align is strongly 
determined by their mutual interactions which provide the key to understanding the 
emergent behaviour at high particle density. In this regime, the interplay between 
microscopic self-motility and anisotropic volume-exclusion interactions leads to complex 
spatio-temporal behaviour that can be directly visualised in two spatial dimensions, i.e. 
for particles moving in planar confinement. 

Quasi two-dimensional systems of self-propelled particles can be realised in a 
number of ways. Autonomously navigating bacteria and other microbes can be confined 
to free-standing thin films [9], between solid surfaces [12] or a liquid-gas interface 
[8J [13] . On larger length scales, active systems can be realised by polar granular rods 
on a flat vibrating surface [TH ITS] or pedestrians moving in complex environments 
[15] . Last not least, colloidal dispersions constitute ideal model systems not only for 
investigating passive matter [TTJ US] but also for active matter composed of self-motile 
colloidal particles. Over the past decade, a number of distinctly different realisations 
of active colloidal particles have been proposed. These include Janus particles driven 
by catalytic processes [HJ [20] or thermophoretic [21] gradients, particles propelled by 
artificial flagella [22] and surface waves [23j [21] driven in an external magnetic field. 
Rather than being spherical most of these particles have an anisotropic rod-like shape 
which is found to play a crucial role in determining the spatio-temporal behaviour of 
active particles [25, 26J. Confining systems to quasi-planar geometries allows for a direct 
visualisation of the particles by means of real-space microscopy and provides fascinating 
opportunities to study the single-particle and collective behaviour of micro-swimmers. 

In this paper we use computer simulation to study a simple model for suspensions 
of rigid, self-propelled rods (SPR) that interact via a Yukawa-segment potential [27, 28J. 
The potential allows for a realistic description of the strong mutual short-range repulsion 
that prevents particles from overlapping. Self-motility is imposed by introducing a 
constant propulsion force along the main orientation axis of each rod. Consequently, 
when two neighbouring active rods collide they align and the aligning force plays an 
essential role in the formation of flocks of coherently moving particles [29]. In our 
study we focus on the collective behaviour of dense suspensions of strongly interacting 
particles and characterise the emergent states by analysing different correlation functions 
as dynamical diagnostics. In order to retain a generic framework we consider the 
overdamped frictional dynamics of the many-body system where the equations of motion 
arise from a simple force balance between the Stokesian frictional force, the collision force 
and the active force on each rod. Likewise, the rod orientations propagate via a torque 
balance involving the frictional and interaction torque acting on each particle. Other 
forces due to e.g. many-body hydrodynamic interactions or thermal fluctuations exerted 
by the embedding solvent are neglected. This allows us to simplify the microscopic 
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equations of motion in such a way that the rod aspect ratio and density constitute the 
main variational parameters of the model. The microscopic self-propulsion force can 
be appropriately scaled out and subsumed into an (effective) Yukawa amplitude which 
only has a weak impact on the emergent behaviour. 

Despite its simplicity the model is capable of predicting a wealth of different steady 
states that hitherto could not be realised within a single framework. Amongst the 
various states we identify an incoherent, disordered dynamical phase at small particle 
aspect ratio and a cooperative swarming state at larger rod anisometry as found in a 
number of particle-resolved models [291 [25j |28j [30]. At high densities and small aspect 
ratios, we find a jammed phase with distinct local crystalline order. This state is rather 
common for passive systems [31] but less obvious for active systems. At large aspect 
ratio and high density, stratified patterns emerge consisting of lanes driven in opposite 
directions. These structures are reminiscent of laning patterns observed for mixtures of 
passive particles (i.e. with no internal driving force) driven in a macroscopic external 
field [32} [331 EH EH]. A similar phenomenon was unveiled recently in mixtures of active 
and passive rod-like particles [36] . For intermediate densities and aspect ratios, we 
find distinct chaotic states characterised by meso-scale turbulent flow patterns with a 
significant vorticity in the velocity field [37]. This type of active turbulence has been 
observed in microbial suspensions [HI EH EH] • Contrary to traditional turbulent flow 
observed at high-Reynolds-number passive fluids the vortices that make up the turbulent 
flow patterns have a uniform mesoscopic size irrespective of the density or particle shape. 

In principle, the full variety of different emergent states advanced here should 
be verifiable for bacterial systems and artificial rod-like colloidal or granular micro- 
swimmers. In a recent study, the statistical properties of the turbulent states as 
predicted from the SPR model have been systematically compared with flow-field data 
of confined bacterial systems [31]. It would be interesting to pursue a more systematic 
comparison with bacterial systems and assemblies of man-made micro-swimmers in order 
to verify the full topology of the predicted phase diagram. 

The remainder of this paper is organised as follows: in section 2 we specify our 
model for self-propelled rods, the corresponding equations of motion and the simulation 
methodology. Numerical results on the non-equilibrium phase diagram are presented 
and analysed in section 3. We conclude in section 4 with a brief discussion of 
possible extensions of the model and we highlight opportunities to observe the predicted 
behaviour in experiment. 

2. Frictional dynamics of a self-propelled-rod (SPR) model 

One of the simplest ways to envisage a suspension of active mesogens is by considering 
a collection of rigid, self-propelled rods each moving with a constant self-motile force F 
directed along the main rod axis (see figure [I]). Mutual rod repulsion is implemented by 
discretising each rod into n spherical segments and imposing a repulsive Yukawa force 
with characteristic decay length A between the segments of any two rods, such that A 
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Figure 1. Coarse-grained representation of a pair of rod- like micro-swimmers with 
n — 5 repulsive Yukawa segments and aspect ratio a = Self-propulsion is provided 
by a constant force F acting along the main rod axis indicated by the orientational unit 
vector u. The total rod pair potential is obtained by a sum over all Yukawa segment 
pairs with distance and is a function of the centre-of-mass distance vector Ar and 
orientations (equation 0). 



defines the effective diameter of the rod of length £ [27J. If two sufficiently long rods 
perform a pair-collision, the interaction results in an effective nematic (apolar) alignment 
while the centres-of-mass attain a certain minimal distance due to the repulsive Yukawa 
forces. The potential energy of a rod-pair a and /3 with orientation unit vectors {u a , u^} 
and center-of-mass distance Ar aj g, is given by 



the distance between the ith segment of rod a and the jth segment of rod (3, with 
h G [—(£ — A)/2, {t — A)/2] denoting the position of segment i along the symmetry axis 
of the rod a. The screening length A defines the effective diameter of the segments such 
that we may introduce an aspect ratio a = ij\ to quantify the effective anisometry of 
the SPR. The case a = 1 corresponds to a single Yukawa point particle {n — 1). For 
a > 1, the number of segments per rod is fixed as n = 3 for 1 < a < 3 and n = \_9a/8] 
for a > 3 with [_•] denoting the nearest integer. 

We focus on the dynamical regime relevant to micro-organisms and artificial self- 
motile colloidal mesogens and we assume the motion of the SPRs to be overdamped due 
to solvent friction (in the zero Reynolds number limit). Since we are interested in the 
collision-dominated dynamics in dense suspensions, we disregard thermal and intrinsic 
fluctuations of e.g. bacterial orientation Consequently, the equations of motion for 
the center-of-mass r a (t) and orientation u a (t) of each SPR are entirely deterministic 




(1) 



where Uq is the potential amplitude, A the screening length, and 
rff = |Ar Q/3 + (liU a - ljUp)\, 



(2) 
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and can be written compactly as 

f T -d t r a = -V ra U + Fu a , (3) 

f R -d t u a = -Vu a U. (4) 

Here, F is a constant self-motility force acting along the longitudinal axis of each rod 
(figure [TJ), U = (1/2) J2p a-p^a Ua/3 the total potential energy, Vu denotes the gradient 
on the unit circle, and 

fr = fo [/||U a U Q + - u a u a )] , (5) 

f* = /o/Rl, (6) 

are the translational and rotational friction tensors (I is the 2D unit tensor) with 
a Stokesian friction coefficient / . The dimensionless geometric factors {f\\, f±, fn} 
depend solely on the aspect ratio a, and we adopt the standard expressions for rod- 
like macromolecules, as given in Ref. jlQ] 

2 " = In a- 0.207 + 0.980a" 1 -0.133a" 2 , (7) 



/ll 

47T 

ira 2 



In a + 0.839 + 0.185a" 1 + 0.233a" 2 , (8) 

In a - 0.662 + 0.917a" 1 - 0.050a" 2 . (9) 

It is expedient to multiply equation ([3J with the inverse matrix f^ 1 : 

d t r a = v u a - f^ 1 • V Ta U, (10) 

where 

^0 = -f^r, (11) 
JO J || 

defines the self-propulsion velocity of a non-interacting SPR. 

In our simulations, we have adopted characteristic units such that X = 1, F = 1, 
and /o = 1, which means that distance is measured in units of A, velocity in units 
of F/ fo, time in units of r = Xfo/F and energy in units of FA. Upon rescaling to 
dimensionless coordinates, three relevant system parameters remain: The dimensionless 
Yukawa amplitude Uq = Uq/ {FX), which determines the hardness of the rod interactions 
relative to their characteristic propulsion energy, the aspect ratio a, and the effective 
volume fraction of the system 



N 
A 



X(£-X) + 



ttA 



(12) 



where the term between brackets denotes the 2D volume A ro( j of a spherocylindrical rod. 
For steeply repulsive Yukawa interactions, the general dynamical behaviour resembles 
that of hard rods and does only weakly depend on the Yukawa amplitude, and we fix 
Uq = 250. The remaining quantities, the rod shape a and volume fraction constitute 
the main steering parameters for our investigations. We simulate the evolution of the 
many-body SPR model as a function of time r = t/ro in a square box of length L 




Figure 2. Schematic non-equilibrium phase diagram of the 2D SPR model at variable 
aspect ratio a and effective filling fraction <j). Values exceeding unity are, in principle, 
possible due to the softness of the Yukawa interactions. The area relevant to self-motile 
bacteria is highlighted in red. A number of distinctly different dynamical states are 
discernible as indicated by the coarse-grained maps of the velocity field v(r,t) (upper 
panels) at time t and the corresponding scalar vorticity field w(r, t) = [V x v(r,i)] • e z 
(lower panels) expressed in units of Tq 1 . 



with periodic boundary conditions at volume fractions in the range 0.05 < (ft < 0.9. 
The simulations are carried out using a time discretisation At = 0.002p~ 1 / 2 , where 
p = NX 2 /A with typically N = 10 4 and rods per simulation. Initial configurations, 
generated from a rectangular lattice of aligned rods with u pointing randomly up or 
down are allowed to relax during an interval r = 1000 before statistics is being gathered 
over an interval r = 20L with L = (iV/p) 1 / 2 the dimension of the simulation box (in 
units of A). Velocity vector fields v(r, t) are constructed by measuring the average 
centre-of-mass velocity within small sub-cells centered around the position r. To this 
end we project the particle positions onto a 2D cubic grid | 1 < i, j < G} and 

measure the average velocity v(t;i,j) in each bin at a given time t. In order to 
test for finite size effects, we consider two different system sizes: 'small' systems with 
N — 1 • 10 4 particles and 'large' systems with N = A ■ 10 4 particles at the same filling 
fraction (ft. The coarse-graining parameter G is chosen adaptively such as to ensure 
each bin to represent the average velocity of about 10 SPRs. Generally, we observe that 
the dynamical structure and order parameters of the emergent states are robust with 
respect to changes in the particle number N, provided N is at least of (9(10 4 ). 



3. Results 



3.1. N on- equilibrium phase diagram for the SPR model 

Upon varying the effective volume filling fraction (ft and the rod aspect ratio a a 
number of qualitatively different dynamical phases emerge. A schematic non-equilibrium 
phase diagram, shown in figure [2j illustrates the importance of the SPR anisometry in 
determining the stationary dynamical state of the system. The low density regime is 
generally characterised by disordered dynamics with little or no cooperative motion. 
Beyond a certain threshold density cooperative motion becomes manifest and translates 
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Figure 3. Evolution of the average SPR velocity as a function of filling fraction <j> for 
a number of particle aspect ratios. Shown is the average velocity component Uii along 
the main rod orientation (a) and the average perpendicular component v± (b), both 
expressed in units of the velocity vq of a free SPR. 



into dynamical states whose structure depends on the intrinsic 'aligning force' of 
the SPRs. Short rods generally jam at high packing fractions whilst very long rods 
(a > 13) exhibit swarming behaviour with large spatio-temporal density fluctuations. 
The swarming and laning phases adjoin a large region of bio-nematic and turbulent flow 
characterised by vortices and extended nematic jet-like structures [37J HI]- 

Generally, the transitions from the dilute phase to regimes with strong cooperative 
motion can be localised by the 2D Onsager overlap density [42J, defined as the density 
corresponding to a singe rod occupying an average area equal to its excluded area 
A ex = (2/tc)(£ — A) 2 + (7r/4)A 2 . The latter expression can be derived from the rod 
dimensions in figure [T] by assuming a pair of spherocylindrical rods with isotropic 
orientations. By combining terms one arrives at the following expression for the overlap 
density: 

,* = Aod = l + 4(a- 1)/tt 

9 A ex l + 8(a- 1) 2 /tt 2 ' 1 ' 

This quantity delimits the density regime beyond which many-body rod collisions 
(exceeding the pair level) are starting to become important and various non-trivial 
emergent states arise. In the sections below we shall discuss these in more detail. 



3.2. Short rods: active jamming 

For small aspect ratios a > 3 a distinct transition towards a jammed state is observed 
upon increasing density. This behaviour is hinted at by the average SPR velocity for 
which we may probe both parallel and transverse contributions via 
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Figure 4. Mean-square displacement of the centre-of-mass for SPRs with aspect ratio 
a = 3. The snapshots depict velocity fields (upper panels) and the SPR coordinates 
(lower panels) for two different bulk filling fractions corresponding to a dynamically 
disordered fluid at cj) = 0.70 and a jammed state at (f> = 0.97. Colour coding is used to 
indicate the orientation u y = u • e y of each rod. 



V± = N ( S P " UaUatt)) ■ Va{t)j , (14) 

where the brackets (■ ■ ■) denote a time average. The results are depicted in figure [3] 
In general, the average parallel velocity decreases monotonically with density as the 
particles get progressively hindered in their motion due to mutual rod collisions. For 
small a the mobility drops rapidly for larger until reaching a threshold level at virtually 
v\\ ~ indicating dynamical arrest. This behaviour is more clearly reflected in the mean- 
square displacement (figure [4]) where a sharp drop in the mobility (over nearly 2 orders of 
magnitude) at <fi = 0.84 marks the onset of jamming. Throughout the density range the 
motion is observed to be sub-ballistic at long times with (r 2 ) ~ r i-75±o.i j am ming 
point depends strongly on particle anisometry as indicated in figure [2] with a marked 
shift towards higher volume fractions upon increasing a. From a structural point of 
view the jamming transition is accompanied by a crossover towards orient at ionally and 
positionally ordered structures as evident from the marked degree of local crystalline 
order at large filling fractions. The velocity maps reveal small pockets of locally enhanced 
particle mobility that point to the presence of dynamical heterogeneities as commonly 
found in glassy systems at finite temperature. A detailed account of freezing and glassy 
behaviour of self-motile spherical Yukawa particles has been reported in Ref. |43j . 
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Figure 5. (a) Enstrophy Q (in units t^ 2 ) versus filling fraction for a number of aspect 
ratios a in the turbulent regime. The maxima correspond to the densities where mixing 
due to vortical motion is the most efficient, (b) Spatial velocity autocorrelation function 
for a number of bulk volume fractions in the turbulent flow regime for two different 
aspect ratios a. 
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Figure 6. Maps of the vorticity field w(r,<) = [V x v(r, t)\ ■ e z expressed in units 
of Tq 1 showing large-scale turbulent flow for SPRs at intermediate aspect ratios: (a) 
(j) = 0.72, a = 5 and (b) = 0.90, a = 7. The snapshots are based on N = 4.10 4 SPRs. 
The lateral box dimensions are 103£ (a) and 78£ (b). 



3.3. Intermediate aspect ratio: vortical states and turbulence 

The maximum in the transverse SPR velocities depicted in figure [3]d suggest that the 
SPRs exhibit some degree of collective swirling motion at moderate densities even at 
small aspect ratios. This type of motion becomes much more manifest at larger a where 
distinct vortical patterns arise akin to turbulent flow. The kinetic energy associated 
with local vortical motion can be measured from the enstrophy per unit area [HJ H3 H6] 
which is defined as: 

n = I (MmTF) , (is) 

where the overbar denotes a spatial average. For slender rods (a > 3) the mean 
enstrophy exhibits a pronounced maximum as a function of the volume fraction <p (figure 
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Figure 7. Power spectra of the kinetic energy for turbulent flow of SPRs with a = 7 
(hi = 2tt/£). Universal scaling behaviour (with scaling exponent —5/3) is observed in 
the intermediate range of wavenumbers k. 



|5p). This maximum signals the density at which vortical motion is maximal. In a 
bacterial suspension this extremum would correspond to the optimal concentration for 
fluid mixing. The range of aspect ratios over which turbulence flow is stable corresponds 
well with the typical aspect ratios of bacterial cell bodies, e.g. a ~ 3 for E. Coli and 
a ~ 6 for B. subtilis (cf. figure [5]). 

The typical size of the vortices that make up the turbulent flow patterns can 
be extracted from the equal-time velocity autocorrelation function (VACF) g v (r) = 
(v(0,t) ■ v(r, £)). This quantity can be obtained from the microscopic SPR coordinates 
{r a , v a } via: 

(E a ~ ka - r^|)(v a • - (v) 2 )) 
9v(r) = \ -^A (16) 

The decay of the VACFs in figure |5]d reveals a typical vortex size of about ~ 5£, 
an estimate that seems rather insensitive to the bulk density and aspect ratio. 
Monotonically decreasing velocity correlations correspond to bio-nematic-type states 
where large-scale nematic jets and vortices coexist |37j whereas negative correlations 
(cf. the curves for a = 7 and > 0.8) represent more pronounced vortical collective 
motion reminiscent of fully developed meso-scale turbulent flow [39]. Typical vorticity 
snapshots are shown in figure |6j 

In order to make a connection with classical 2D turbulence in high- Reynolds number 
fluids we have calculated the energy spectrum which can be obtained as a Fourier 
transform of the VACF: 



E{k) = — I drexp[-ik-r](v(0,t) • v(r,t)). 

2-7T J 



(17) 
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Figure 8. Spatial auto-correlation functions of the vertical velocity component g v (r) 
(equation (16)) at various filling fractions <j> for SPRs with a = 16. Particle snapshots 
of the stationary states at (a) <f> — 0.05 - disordered, incoherent motion (b) <f> = 0.21 
- swarming, (c) 4> — 0.32 - initial stage of laning formation, and (d) 4> — 0.42 - fully 
developed laning state. The lateral box dimension corresponds to 38£. Colour coding 
is used to indicate the orientation u y = u • e y of each rod. 



An alternative and more formal definition reads (v 2 ) = 2 J °° dkE{k) where E(k) reflects 
the accumulation of kinetic energy over different length scales. The results in figure [7] 
suggest asymptotic power law scaling regimes for intermediate fc-values with a power- 
law exponent close to the characteristic fc _5,/3 -decay predicted by the Kolmogorov- 
Kraichnan scaling theory jUl HB] for (passive) 2D turbulence in the inertial regime. In 
the present case, however, inertia is absent on the particle scale because the SPR motion 
is completely overdamped but it is possible that the self-propulsion establishes 'effective' 
collective inertial effects on larger scales which could explain the observed k~ 5 ^ 3 decay. 
Contrary to regular turbulent flow where energy is injected on the macroscopic scale, 
active turbulence is characterised by forcing on the microscopic scale. In general, the 
transport of kinetic energy towards smaller k becomes significantly damped on larger 
length scales [13] as highlighted by the low- A; plateau in the power spectra in figure [7j 
We refer the reader to Ref. [39] for a more detailed discussion comparing meso-scale 
turbulence in active suspensions and regular high-Reynolds-number turbulent flow. 

3.4- Long rods: swarming and lane formation 

At low to moderate density slender rods with a > 10 tend to form large compact 
flocks reminiscent of the cooperative motion observed in large groups of organisms, e.g. 
fish schools, bird flocks (see figure [8^i-b) [1HJ ED]- At larger volume fractions distinct 
laning patterns emerge that consist of cooperative stratified motion (figure [8]>d) . The 
transition from flocking to laning as a function of the volume fraction can be localised 
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from the equal-time VACF for the velocity components along the lane directions (in this 
case the vertical component). The characteristic decay of the VACF g v (r), shown 
in figure [8j allows a distinction between the disordered state at low densities where rod 
clusters are small and velocity correlations decay rapidly and the emergence of large 
cooperative flocks with velocity correlations spanning several dozens of rod lengths. 
A marked divergence of the correlation length occurs around <fi ~ 0.4 where the flocks 
start to span the entire system and self-organise into lanes moving in opposite directions. 
The laned patterns remain stable throughout the sampled time interval and no sign of 
break-up is observed over time even for large systems. We have verified the stability 
of the lanes against small thermal fluctuations of the rod orientations that could be 
induced by the embedding medium or by some internal source, e.g., bacterial flagella. 
The rotational fluctuations are represented by a Gaussian white noise contribution Au 
to the equation of rotational motion of each rod a (cf. equation Q): 

<9 t u a = -f^-V aa [/ + Au a . (18) 

The stochastic term has zero mean (Aii ia ) = and correlations {Au ia {t)Aujp{t')) = 
2D* R 5ij5 a p5(t — t') (with i — x,y) in terms of some effective rotational diffusion rate 
D* R . Although 'run-and-tumble' motion as commonly observed in bacterial systems 
(notably E. Coli (HE!]) is strictly non-Brownian at short times, its long-time behaviour 
is well-captured by a rotational diffusion process with a diffusion constant much larger 
than the Stokes-Einstein value D R = k R T/f f R (where k R T is the thermal energy) 
for passive Brownian rods |51j . The strength of the tumbling motion is conveniently 
expressed in terms of the dimensionless tumbling parameter iD* R j v which is the ratio of 
the translation time a free SPR needs to swim over a distance i and the typical tumbling 
time \ jD* R . Typical values for E. Coli and other swimming bacteria are £D* r /vq ~ 0.01 
[I]. In the dense regime, the particle velocities are dominated by rod-rod collisions 
rather than thermal fluctuations and the intrinsic rotational diffusivity of the SPRs 
does not incur any qualitative change to the laning structures. In general, we assume 
that the spatio-temporal states and the topology of the phase diagram is robust against 
weak fluctuations in the swimming direction of the SPRs. We remark that similar 
laning states were encountered at finite temperature in binary mixtures of SPRs with 
different self- motility [36] . In both classes of driven systems laning instabilities occur 
if the disparity between the species mobility exceeds a certain threshold. In our case, 
however, such an intrinsic driving force is absent since all particle have equal mobility. 

4. Conclusions 

We have studied the collective dynamical behaviour of a simple two-dimensional model 
of self-propelled rigid rods (SPR) by means of numerical simulation. Depending on 
the rod shape and density, the SPR model exhibits a wealth of different emergent 
dynamical states including swarming, turbulence, laning and jamming. Although many 
of these states have been encountered in various set-ups, most notably (mixtures of) 
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spherical particles in different external fields, the SPR model is able to generate these 
dynamical states upon variation of only two basic system parameters; the particle shape 
and density. The present approach may therefore serve as a benchmark to characterise 
the collective properties of different classes of self-motile organisms and artificial micro- 
swimmers of various shapes. As for the turbulent state, it was recently shown that the 
SPR model is capable of reproducing the velocity statistics obtained from experiments 
on strongly confined Bacillus subtilis suspensions [SH]- Future experiments on dense 
systems of self-propelled particles with low and high particle anisometry will hopefully 
allow for similar comparisons for the jammed and laned state, respectively. 

Future efforts could be aimed at extending the SPR model and the associated 
equations of motion by accounting for e.g. multi-body hydrodynamic interactions 
mediated by the solvent, particle flexibility and body forces transmitted by chemical 
gradients (chemotaxis) that could be relevant in concentrated bacterial systems. The 
influence of thermal fluctuations could be incorporated if one wishes to assess the 
effect of translation and rotational noise (bacterial tumbling) in more detail. It is 
also desirable to explore the SPR model in three spatial dimensions, for instance, to 
study the phenomenology of fully developed 3D meso-scale active turbulent flow which 
has been unexplored so far. Finally, it would be challenging to construct microscopic 
theories that are capable of predicting the observed emergent states advanced in this 
study. Dynamical density functional theory for anisotropic particles [2H1 E21 E3] could 
provide a promising avenue for this. 
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